Modelling galaxy halos using dark matter with pressure 
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We investigate whether a dark matter with substantial amounts of pressure, comparable in magni- 
tude to the energy density, could be a viable candidate for the constituent of dark matter halos. We 
find that galaxy halos models, consistent with observations of flat rotation curves, are possible for a 
variety of equations of state with anisotropic pressures. It turns out that the gravitational bending 
of light rays passing through such halos is very sensitive to the pressure. We propose that combined 
observations of rotation curves and gravitational lensing can be used to determine the equation of 
state of the dark matter. Alternatively, if the equation of state is known from other observations, 
rotation curves and gravitational lensing can together be used to test General Relativity on galaxy 
scales. 

PACS numbers: 
Keywords: 



I. INTRODUCTION 

The rotation curves of spiral galaxies are one of the 
most direct probes of dark matter 0- In these galax- 
ies neutral hydrogen (HI) clouds are observed at large 
distances from the center, much beyond the extent of 
the luminous matter. The frequency shifts in the 21 cm 
HI emission from these clouds allow their velocities to 
be measured. These clouds move in circular orbits with 
velocity Vc{r). The orbits are maintained by a balance 
of the centrifugal acceleration jr and the gravitational 
pull GM(r)lr^ of the total mass M(r) contained within 
the orbit . This allows us to use the observed rotation 
curves (wc as a function of r) to determine the mass pro- 
file M (r) = Tv\ (r) /G of the galaxy. 

The rotational velocities are observed to increase near 
the center of the galaxy and then remain nearly constant 
at a value Vc ~ 200km/s. This leads to a mass profile 
M(r) = rv'^/G where the mass within a distance r from 
the center of the galaxy increases linearly with r even at 
large distances where there is very little luminuos mat- 
ter seen. This is usually explained by postulating the 
existence of some dark (invisible) matter distributed in 
a spherical halo around the galaxy. A modified law of 
javity has been proposed as an alternative explanation 
, we do not consider this here. 
There is at present no clear picture as to what con- 
stitutes the dark matter in galaxy halos. The prevalent 
belief is that it is a pressureless medium (Cold Dark Mat- 
ter) possibly made of weakly interacting massive particles 
Q (WIMPS). Here we consider the possibility that the 
dark matter in galactic halos has substantial pressure, 
comparable in magnitude to the energy density. 

The Cold Dark Matter (CDM) model has been ex- 
tremely successful in explaining the observed CMBR 
anisotropics and large scale structures in the universe 
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The evolution of cosmological perturbations is dras- 
tically modified in models where the dark matter is an 
ideal fluid with significant amounts of pressure, and such 
models are ruled out by these observations. Our study 
considers the possibility that the dark matter could be 
a medium which is not an ideal fluid and has significant 
pressure. Such a situation is generic if the dark matter 
is modelled as a scalar field (for references on the use 
of scalar fields in modelling dark matter/dark energy see 
the papers cited in T]), a scenario which has come up 
quite often in the various attempts to solve the prob- 
lems being faced by the CDM model on small scales 6] . 
This, and the fact that there has, till date, been no direct 
detection of the CDM particles motivates us to explore 
the possibility that the dark matter in galaxy halos may 
have pressure. In this work we do not restrict ourselves 
to any particular model for the dark matter. We consider 
a general situation where the pressure inside the halo is 
anisotropic and is related to the energy density through a 
linear equation of state. Other attempts at interpreting 
the flat rotation curves in the recent past are available in 
i- 

We investigate whether it is possible to explain the flat 
rotation curves seen in the outer parts of spiral galaxies 
using dark matter which has pressure. The Newtonian 
analysis of the rotation curves outlined in the beginning 
of this section fails when the pressure is comparable to the 
energy density, and a relativisitic treatment is required. 
In the Newtonian analysis the gravitational field inside 
the halo is described by a single function, the gravita- 
tional potential, which is completely determined by the 
observations of the rotation curve. In the relativisitic 
treatment the gravitational field inside a static, spheri- 
cally symmetric dark matter halo is described by the two 
unknown metric coefficients gao{r) and grrif)- In section 
2 of this paper we investigate to what extent the metric 
inside the halo can be determined by observations of the 
rotation curve. We show that while 500 (j') is completely 
determined by just the rotation curve, inputs regarding 
ernetilel theory of gravity and the dark matter equation of 
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state are required to determine grr{r). 

In section 3. we calculate the gravitational lensing 
properties of the dark matter halo. The usual analy- 
sis of gravitational lensing does not take into account 
pressure. Here we analyse photon geodesies inside the 
halo and use these to calculate the gravitational lensing 
deflection angles. Subsequently, we study to what extent 
the lensing properties of dark matter halos change when 
pressure is introduced. We investigate whether it is pos- 
sible to use combined observations of rotation curves and 
gravitational lensing to constrain the equation of state of 
the dark matter in galaxy halos. In section 4. we discuss 
our results and present conclusions. 



II. THE GEOMETRY WITHIN THE HALO 

It is assumed that spiral galaxies are embedded in 
static, spherically symmetric halos of radius R. The 
gravitational field outside the halo is described by the 
exterior Schwarzschild metric with M = Rv'^/G. The 
gravitational field inside the halo is represented by two 
unknown functions goo{^) = e^*^""^ and grr{^) = e^^^''^ 
which parametrise the metric, and the proper time inter- 
val dr is 
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We first investigate to what extent the gravitational field 
inside the halo can be determined from observations of 
the rotation curves. 

The HI clouds from which the frequency shifted 21 cm 
emission is observed are treated as test particles mov- 
ing in circular orbits, in the plane of the disk of the 
galaxy under the gravitational influence of the halo. It 
is assumed that the HI clouds do not contribute signif- 
icantly to the halo matter. We choose coordinates so 
that the disk of the galaxy is aligned with the equato- 
rial plane of the halo {6 = tt/2). The orbits of the HI 
clouds have two conserved quantities, E — e^'^cdt/dr 
and L = r'^dcp/dr which, for circular orbits, are related 
as = e^* [c2 + L^/r^] . The HI clouds at different 
radii will have different values of E and L, and we de- 
note this by E{r) and L{r). 

A distant observer in the plane of the disk receives HI 
emission from a cloud located along a line of sight which 
is tangential to the orbit of the cloud. The observed shift 
in the frequency of the radiation, when the cloud is mov- 
ing away, is Av « —VeL{r)/{cr)^ provided e^*^'"' ^ 1 
and L{r)/[cr) <C 1. Later in this paper we show that 
these two approximations are valid for our analysis. An 
extra angular dependence comes into the frequency shift 
if the line of sight is at an angle to the direction of motion 
of the HI cloud. Observations of the frequency shifts in 
the HI radiation show that, barring the central parts of 
galaxies, the ratio L{r)/{cr) is nearly constant at a value 
L{r)/{cr) — Vc/c ^ 7 x 10^"*. The ratio Vc/c being ex- 
tremely small, in our analysis, we retain only the leading 
order terms in Vc/c, and ignore all the higher terms. 



We next shift our attention to the equation of motion 
of the HI clouds. For circular orbits the geodesic equation 
reduces to 
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where $ = d^/dr. This can be integrated to obtain 
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This equation allows us to determine the metric coeffi- 
cient .goo('') inside the halo using just the observed rota- 
tion curves. For a flat rotation curve this gives us 
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which gives the metric coefficient goo(r) to be 
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Here the constant of integration in equation (|3Jl has been 
chosen so that the metric coefficient 500 (^) matches the 
exterior Schwarzschild metric at r — R. 

The rotation curve determines only one, namely <I>(r) 
of the two unknown functions $(r) and A(r), which are 
required to describe the gravitational field inside the halo. 
We also see that the metric coefficient goo{r) varies very 
slowly with r, and e^**^''-' ~ 1, the change being less than 
0.5% when r changes from O.Oli? to R. The fact that 
goo can indeed be obtained from the rotation curves, has 
been partially noted in the papers in |^ and in some of 
the references on halo dark matter 0. 

Next we address the question ~ how to determine A(r)? 
We proceed by solving Einstein equations for gravity 
Gfj-u = {8nG/c'^)T^„ to determine A(r) inside the halo. 
Solving this requires the equation of state (or equiva- 
lently the energy momentum tensor T^j/) for the matter 
which makes up the halo. We assume that the dark mat- 
ter which makes up the halo is a fluid with energy den- 
sity c^p(r), radial pressure Pr{r), and tangential pressure 
Prii") in the 9 and <j) directions. Also, the gravitational 
field inside the halo is weak ie.^{r) ~ A(r) ~ {vc/c)"^ 
which allows us to linearise the Einstein equations retain- 
ing only terms linear in (wc/c)^. Inside the halo, using 
equation ^ for <i>(r), the Einstein equations reduce to 
three equations for A(r) 
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which have to be solved with the boundary condition 
A(i?) = (vc/c)'^ and the requirement p{r) > 0. 



3 



In the absence of any pressure {Pr 
recover the usual Newtonian solution 
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This solution corresponds to the singular isothermal 
sphere which produces a flat rotation curve, We have 
solved equations ©, {Tj) and © under various assump- 
tions on the equation of state. It is convenient to express 
the results in terms of Aat and pN (r) which arc the soul- 
tions in the absence of pressure. 

The most obvious possibility is to consider a medium 
with isotropic pressure {ie.Pr = Pt)- We find that the 
only solution where the metric matches the Schwarzschild 
metric at r = i? and the energy density is positive has 
P^ = Prp = 0. 

We next consider two possible anisotropic equations 
of state A. Pr = Wrpc^ where we use eq. © to deter- 
mine Pt, and B. Pt — wtPC^ where we use eq. (O 
to determine Pr- As mentioned in the Introduction, 
such anisotropic pressures can indeed arise in the energy- 
momentum tensor of a real or complex scalar field. We 
recall that the energy momentum tensor for a real scalar 
field in a potential V{4)) is given as : 



When the gradient terms are small, and the scalar field 
is frozen at a potential minimum (pmin 7^ 0, it generates 
isotropic, negative pressures P = —pc^ = —V{(j)min)- On 
the other hand, the pressure is anisotropic i.e. the ra- 
dial pressure Trr is different from the tangential pressure 
Toe, if the radial derivative of the scalar field makes a 
significant contribution to the energy momentum tensor. 
Various scalar field configurations, including some with 
anisotropic pressures (for example, Matos and Nunez, 
astro-ph/0303455 in 7]), have been used extensively in 
cosmology to model the dark matter and the dark energy. 

We now discuss the results for the two cases (A) and 
(B) mentioned above. 



A. /V = WrpC 

We find that solutions exist for Wr'>—1. If we match 
match A(r) to the Schwarzschild metric at the outer edge 
of the halo, the density turns out to be negative for Wr < 
— 1 and this range is ruled out. The allowed solutions are 
tabulated below. 



Tf^ = d,cj)d,(b ~ -.9,, [d'^^bdk^ + 2V{(b)] (10) 
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l+w, 



1 + Wr \B. 



1 + 1(77- 



In 



Prir) = (cV2)p^(r)x 
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B. Pt = WTpc? the density is negative if wt < —1/2. The solution for 

the allowed range of wt^s shown below 

We find that solutions exist for wt > —1/2. For wt = 
— 1/2 we have an absurd situation where A(r) = 0, and 
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III. GRAVITATIONAL LENSING 

We now move on to calculate the bending of light by 
the halos. A light ray which goes past the halo without 
entering it propagates entirely in a Schwarzschild metric. 
The light ray is deflected by an angle A = AGM/c^h = 
A{vc/ cYiR/b) where b is the impact parameter, and in 



this case b > R. The bending of a light ray which passes 
through the halo is determined by the metric inside the 
halo and this depends on the equation of state of the dark 
matter. 

The light ray is assumed to lie in the = 7t/2 plane. 
NuU-geodesics in the halo metric are completely deter- 
mined by the two conserved quantities E and L defined 



4 



in section 2. The trajectory of the light ray is determined 



by the following equation QJ^ for «(</>) — 



du 
14, 



^ g-2(<I.+A) 



(11) 



Retaining only terms to order ^ (v /cf' , and using 1/6^ 
E"^ /L"^ we have 
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It is convenient to replace w by a new variable y where 
y = u + a{y) and 



a(y) 



($ + A) /&2 - Ay2 



We can now integrate H12() to obtain 
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In equation p4l) the limit y = is when the in-coming 
photon is very far away from the center of the halo, y — 
1/b is when the photon is closest to the center of the halo 
and (j) is the change in the angle subtended by the photon 
between these two positions. Setting a{y) = allows 
us to calculate the photon trajectory in the absence of 
any gravitation field, whereas the term involving a is the 
change produced due to the gravitational field of the halo. 
This gives the deflection angle to be 



A = -2 



'/"da 
dy 



y 



dy 



(15) 



where the factor of 2 arises from the contributions of the 
in-coming and the outgoing photon trajectory. 

The integral in eq. (|15|l is evaluated in two parts where 
we use the Schwarzschild metric to calculate a outside 
the halo ( l/y > R) and we use the solutions obtained in 
Section 2. inside the halo {l/y < R). The results for the 
deflection angle A are presented in units of [vcjcY using 
a variable 5 = /S./{vc/cY. We recover the familiar result 
5 = 4(i?/6) corresponding to the exterior Schwarzschild 
geometry when the impact parameter is larger than halo 
radius {b > R). In the situation when the light ray does 
enter the halo {b < R), eq. l(T3|l can be evaluated analyt- 
ically using Mathematica for both the equations of state 
considered in Section 2. The resulting expressions are 
rather lengthy and we do not present them here. In the 
absence of pressure we recover the familiar result d — 2Tr 
for an infinite, isothermal halo. Wc recover this result 
from eq. IjlSI) in the limit {b/R) — > 0. This can be in- 
terpreted as either the halo raduis becoming very large 
for a fixed impact parameter, or the impact parameter 
approaching very close to the center of the halo with the 
halo radius remaining constant. The results for the de- 
flection angle in the presence of pressure, under the two 
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FIG. 1: This shows the deflection angle A = {vc/c)^5 as a 
function of the impact parameter b for a light ray passing 
through the dark matter halo of a galaxy which has a flat 
rotation curve with speed Vc- The equation of state of the 
dark matter is assumed to be P,- = Wrp(? and the values of 
Wr for which the deflection angles have been calculated are 
shown in the figure. 
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FIG. 2: This figure is the same as Figure 1 except that equa- 
tion of state of the dark matter is assumed to be Pt = wrpc^ 



different equations of state assumed in section 2. are 
shown in figures 1 and 2. 



IV. RESULTS AND CONCLUSIONS 

We have addressed the question if the dark matter in 
the halos of spiral galaxies could have significant amounts 
of pressure, which in general could be anisotropic. This 
requires a relativistic treatment where the gravitational 
field inside the halo is described by the two metric coeffi- 
cients goo{r) and grr{r). We find that observations of the 
rotation curve fully determine goo(^) irrespective of the 
equation of state of the dark matter. As a matter of fact 
goo (f) is determined without any reference to the theory 
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for gravity, the only assumption being that the gravita- 
tional effects can be represented through the space-time 
metric. 

Determining grr{i') requires us to assume a theory for 
gravity, Einsteins's theory of General Relativity in this 
case, and an equation of state for the radial pressure 
Pr and the tangential pressure Pt of the dark mat- 
ter. We consider two possibihties (A^) Pr = Wrpc^ 
with Pt being determined from Einstein's equations and 
(B.) Pt = WTpc^ with P^ being determined from Ein- 
stein's equations. It may be noted that there is no solu- 
tion matched to the Schwarzschild exterior metric at the 
boundary of the halo (r = R) if the pressure is assumed 
to be isotropic {Pr = Pt)- 

Solutions matched to a Schwarzschild exterior solution 
exist for (A.) Wr > -1 and (B.) wt > -1/2. Mak- 
ing use of the calculated grri^) we determine the density 
and pressure profiles for the halos. The density profile 
is shallower than the Newtonian density profile only 
for -1/2 < WT < 0. 

For Pr = ~pc^ we have an interesting situation where 
<?rr(f) — — 5'oo^(^) both inside and outside the halo. The 
density also goes to zero at the boundary, unlike in all the 
other cases. It should be noted that we do not attempt 
to match the density and pressures at the boundary of 
the halo, and in fact it is not possible to do this for a flat 
rotation curve and a Schwarzschild exterior. 

We next investigate if it is possible to constrain grr 
or equivalently the equation of state using gravitational 
lensing. The bending of light rays was calculated for the 
cases (A.) and (B.) and the results are summarised below. 

(A.) For I |> the bending angle is always larger 
than the case where there is no pressure. When Wr > 0, 
the deflection angle progressively increases as the halo 
radius R is increased keeping the impact parameter b 
fixed, and it diverges in the limit (b/R) 0. On the 
other hand, for Wr < the behaviour of the deflection 
angle is similar to the situation when there is no pressure, 
except that in the limit {b/R) ^ it saturates at a value 



larger than A^r = 2tt{vc/c)'^ which is the deflection for 
an infinite halo with no pressure. 

(B.) The generic features of the deflection angle are 
different for the two cases wt > and wt < 0. For wt > 
0, the deflection angle is larger than Ajv, and it increases 
progressively as the halo radius is increased, diverging in 
the limit {b/R) — > 0. For wt < 0, the deflection angle is 
less than A^r and it converges to around CSA^r in the 
limit {b/R) 0. 

We find that in case (B. ) with wt < 0, the only sit- 
uation where the deflection angle is less than Ajv, the 
density profile is shallower than r^^. The density has a 
r^^ behaviour or is steeper in all the other the cases. It 
may be worthwhile to remind the reader here that in all 
the cases considered here both the energy density and 
the pressure contribute to the bending of light. We next 
ask if there is any pattern seen in the form of grr{r) (or 
effectively A(r)) and the behaviour of the deflection an- 
gle. We find that wt < is the only case where X{r) is 
a monotonically increasing function of r. This leads to 
a picture where the gravitational deflection angle is Ajv 
in the absence of pressure where A(r) is a constant, the 
deflection is more than Ajv if A(r) is a monotonically de- 
creasing function of r and the deflection is less than 
if A(r) is a monotonically increasing function of r. 

In conclusion, the flat rotation curves observed in the 
outer parts of spiral galaxies can be equally well ex- 
plained by a halo which is made up of dark matter with 
anisotropic pressures. The bending of light rays pass- 
ing through the halo is found to be highly sensitive to 
the pressure content of the dark matter. This holds the 
possibility that combined observations of rotation curves 
and gravitational lensing can be used to determine the 
equation of state of the dark matter in galactic halos. Al- 
ternatively, if the equation of state of the dark matter is 
determined from other observations, then combined ob- 
servations of rotation curves and gravitational lensing can 
be used to test Einstein's theory of General Relativity on 
galaxy scales. 
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